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ABSTRACT 

A recently proposed k — e model for low Reynolds number turbulent flows (Yang and 
Shih, 1991) was modified by introducing a new damping function The modified model is 
used to calculate the transitional boundary layer over a flat plate with different freestream 
turbulence levels. It is found that the model could mimic the transitional flow. However, 
the predicted transition is found to be sensitive to the initial conditions. 


1. Introduction 

In a quiescent environment, transition is initiated by the amplification of Tollmien- 
Schlichting waves. These waves eventually breakdown, giving rise to turbulent spots, 
which are precursors of the turbulent boundary layer. While in an environment with high 
freestream turbulence, say the flow passing over a turbine blade, turbulent spots are formed 
due to the transport of turbulence from the freestream to the boundary layer, rather than 
from the T-S wave amplification. 

Priddin (1975) was the first to show that the low Reynolds number two equation mod- 
els have the potential to predict transitional flows under the influence of the freestream 
turbulence. A detailed calculation procedure was given by Rodi and Scheuerer (1985), in 
which the Lam & Bremhorst low Reynolds number k — e model was used. More recently, 
a comparative study of the performance of existing low Reynolds number Jb — e models in 
predicting laminar-turbulent transition was made by Fujisawa (1990). 

Recently, a low Reynolds number k — e model was proposed by the authors (Yang and 
Shih, 1991). The proposed model was found to perform quite well for turbulent channel 
flows and fully developed turbulent boundary layer flows at different Reynolds numbers. 
The purpose of this work is to test the capability of the proposed model in predicting 
the laminar-turbulent transition. The test cases are the transitional boundary layer flows 
over a flat plate with a freestream turbulence level of 3% (Case T3A) and 6% (Case T3B) 
respectively. These cases are the benchmark cases in an ongoing project coordinated by 
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Savill (1991), testing the capability of turbulence models in predicting transitional flows. 

The plan of the paper is as follows: in Section 2, we list the equations to be solved, 
and the corresponding initial and boundary conditions; in Section 3, we brief the numer- 
ical aspects of the calculation; the results of the calculation are shown in Section 4 and 
conclusions are presented in Section 5. 


2. Mathematical Formulations 


The boundary layer approximation is used. Upon using the eddy viscosity assumption, 
the equations for the mean field over a flat plate boundary layer are: 


dx dy 


( 1 ) 


uW + v °Z = ± 

dx* dy dy 


( i / + vt) 


dU 

dy 


( 2 ) 


where x,y are the coordinates along and normal to the plate, and U, V are the mean velocities 
in the x, y directions, respectively. 

Recently, Yang and Shih (1991) proposed a low Reynolds number k — e model to close 
the above equations. In this model, the eddy viscosity is given by 

v T = cJ^kT (3) 

where the time scale T has the Kolmogorov time scale as the lower bound: 

k 


T = - + (-) 1/2 


(4) 


The transport equations for k and e are 


v » + v «i = » 

dx dy dy [ 


v T dk 

[V + )a- 

CT/i, dy 


, dU . , 

+ M-^) -« 
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where 


U— + V— - — 
dx* dy dy 


v T de 

{? + — )q~ 

<r t dy 


+ 


Civt(-q ^) 2 ~ Cifot 


Rl 


+ E 


h = 1.0- 0.22 exp(-j|), 


E 


, dU . 2 
vu T {^) , 


( 6 ) 


and the turbulent Reynolds number is defined by Rt = — 


2 


The model constants used are the same as those in the Standard Model, i.e. 
Cp = 0.09, Ci — 1.45, Cj = 1.92, Uk = 1.0, <r c = 1.3. 


Thus, the performance of the model in high turbulent Reynolds number flows is ensured. 

The damping function / M used in Yang and Shih (1991) is a function of y+. It has been 
documented (see Stephens and Crawford, 1990, for example) that the y + dependence in the 
damping function / p would lead to an early transition in the model prediction, and a lack 
of minimum c/ (which signifies the start of transition) when the freestream turbulence is 
high. We here introduce a new damping function which is a function of Ry, where R v is 
defined by 

(7) 

The form of is chosen to be 

U = [1 “ exp(-aiR v - a 3 Rl - a 6 R *)] 1/2 ( 8 ) 

where the constants in the damping function take the following values 

a.\ = 1.5xl0~ 4 ,a 3 = 5.0xl0 -7 , a$ = l.OxlO -10 . 

These values are determined by calibration with the DNS data for channel flows. The 
damping function chosen gives the correct asymptotic behavior near the wall. Far away 
from the wall, the damping function approaches one as required. The shape of the damping 
function is shown in Fig. 1. 

The above system of parabolic equations need to be supplymented by the boundary 
conditions at the wall and at the freestream, and the initial conditions at the starting point 
of the calculation. At the wall 


U = V = k = 0, 

= i ~dT ) ’ 


(9) 


At the edge of the boundary layer, the flow variables are given by the freestream values, 


i.e. 


U - f7 e} k — fc e ,€ = e c . 


( 10 ) 


The Jb e , e e in the above are found from the transport equations for k and e, with the condition 
that the gradient of the flow variables in the y direction vanishes as the free stream is 
approached. Thus, 

dk e 
dx 

dc t _ C 2 e c 
dx T ' 


U e - 

u e - 


= 


3 


( 11 ) 

( 12 ) 


Jfe e o and e e o (i.€^ the values of k e and e e at the leading edge) are needed. k e q is obtained 
from the experiment , and e e o is determined in such a way that the resulting k e (x) profile 
agrees with the experiment. In the cases considered here, 

^ = 1.18ll0 -3 , JJfJT = 2.86x10“* 

HI U*/L 

for T3A, and 

= 4.94xl0“ 3 , jjjjY = 1.22xl0- 5 

V} U?/L 

for T3B, where L = 1 mm is the reference length scale. 

A major difficulty in calculating the transitional boundary layer through these parabolic 
equations is the lack of information on the initial conditions, particularly, the initial con- 
ditions for k and e. Following Rodi and Scheuerer (1985), we adopt the following initial 
conditions. At the starting location where the boundary layer is still laminar, 

jj- = V B u * («.,* = = ( 13 ) 

with a further constraint that e > e e . 


3. Numerical Aspects 

i 
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A semi-implicit finite difference scheme is used to solve the momentum equation and 
the transport equations for k and e. The coefficients for the convective terms are lagged one 
step in the marching direction, and the source terms in the k and e equations are linearized 
in such a way that numerical stability is ensured. 

A variable grid spacing is used to resolve the sharp gradient near the wall. The grid 
distribution is controlled by Syt/ 6yi—i = ct . Both ot and the total number of grid points, N, 
are varied to ensure the grid independence of the numerical results. The marching step size, 
Sx, is also varied to ensure accuracy. It is found that N = 150, a = 1.03, and Sx = 0.02 are 
sufficient to give a solution with less than a 1.0% error. 

1 4. Results of Calculations 

1 

The modification of Yang and Shih (1991) is' first tested for fully developed turbulent 
flow passing a flat plate boundary layer at Rc$ — 1410, where the DNS data is available 
(Spalart, 1988). In the freestream, the turbulence is assumed to be zero. Arbitrary profiles 
^ for the initial conditions are used. The results for the mean velocity U and the turbulent 

kinetic energy Jfe are shown in Fig. 2 and Fig. 3, respectively. Also shown in these figures 

j 
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are the DNS data. The agreement is very good. Same agreement is also obtained for the 
turbulent boundary layer at Re$ = 670 and the turbulent channel flow at Re r = 395. 

Next, calculations are made for the transitional boundary layer passing a flat plate 
with a freestream turbulence level of 3% (Case T3A) and 6% (case T3B), respectively. 
The calculations are started at Re$ — 100 where the boundary layer is still laminar. The 
distribution of the shape factor H as a function of Re x is shown in Fig. 4. As Re x increases, 
H varies from 2.59 (the value for the laminar boundary layer) to 1.41 (the value for the 
turbulent boundary layer.) The change is more gradual for the high freestream turbulence 
level (Case T3B) than for the low freestream turbulence level (Case T3A), in agreement with 
the experiment. Fig. 5 shows the distribution of the skin friction coefficient c/asa function 
of Re x for Case T3A. Also shown in this figure is the experimental result reported in Savill 
(1991). It is found that Cf varies from a laminar relation to a turbulent relation. Thus, the 
model prediction does mimic transition. However, the predicted onset of transition (defined 
as the point of minimum c/) occurs earlier than that found in the experiment. The model 
prediction gives Re x = 8.0xl0 4 while the experimental result is Re x = 1.3xl0 5 . Similar 
comparison for Case T3B is presented in Fig. 6. In this case, the predicted transition starts 
at Re x = 4.2xl0 4 while transition starts at Re x = 5.9xl0 4 in the experiment. The effect 
of freestream turbulence on transition is captured by the model. It is seen that the higher 
the freestream turbulence, the earlier the transition. In both cases, the predicted transition 
from laminar boundary layer to turbulent boundary is sharper than the experimental one. 

Since it is unlikely that the initial conditions used are the physical profiles for a tran- 
sitional boundary layer at Re$ = 100 found from an experiment, the effect of the initial 
conditions has to be studied before the results of the calculations can be assessed. Since we 
are solving a system of parabolic equations, the change in initial conditions is equivalent to 
a change in the location of the starting point with the same initial profiles. Calculations 
were carried out with different starting locations. Fig. 7 and Fig. 8 show the c/ distribution 
obtained from starting the calculation at Re$ = 10 for Case T3A and Case T3B, respec- 
tively. It is found that the predicted transition hew moved upstream. Transition starts at 
Re x = 3.2xl0 4 for Case T3A and Re x = 2.0xl0 4 for Case T3B. Thus, the initial conditions 
have a significant effect on the results of the calculation. This finding is in agreement with 
that reported in Schmidt and Patankar (1988). 


5. Conclusions 

A recently proposed low Reynolds number Jfc — e model is modified through the intro- 
duction of a new damping function which depends on Ry rather than y+. This modification 
is found to perform well for both fully developed turbulent flat plate boundary layer flows 
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and channel flows. It is also found that it has the capacity to capture the transition from 
a laminar boundary layer to a turbulent boundary layer. However, the predicted transition 
is found to be sensitive to the initial conditions, or equivalently, the starting location of the 
calculation if the same initial profile is used. 

Acknowledgement } 

The DNS data used was kindly made available to us by Dr. N. Mansour of NASA Ames ■ 

Research Center. The first author (ZY) appreciates the help provided by Drs. C. Speziale j 

of ICASE and S. Thangam of Stevens Institute of Technology during the ICASE/LaRC •' 

Workshop on Transition and Turbulence. 

References ■ 

Fujisawa, N. 1990 “Calculation of transitional boundary layers with a refined low Reynolds | 

number version of a J fe - e model of turbulence,” Proc. Engineering Turbulence Modeling | 

and Experiment, pp23-32. | 

Priddin, C.H. 1975 “The behaviour of the turbulent boundary layer on a curved porous | 

wall,” PhD Thesis, Imperial College, London. | 

Rodi, W. and Scheuerer, S. 1985 “Calculation of laminar-turbulence boundary layer tran- I 

sition on turbine blades,” Proc. 65th AGARD-REP Sym. on Heat Transfer and Cooling in g 

Gas Turbine, Bergen, Norway, ppl8-l to 18-12. 

Savill, A.M. 1991 “Turbulence model predictions for transition under free-stream turbu- 
lence,” RAeS Transition and Boundary Layer Conference, Cambridge, England. j 

Schmidt, R.C. and Patankar, S.V. 1988 “Two-equation low-Reynolds number turbulence “ 

modeling of transitional boundary layer flows characteristic of gas turbine blades,” NASA _ 

CR-4145. | 

Spalart, P.R. 1988 “Direct simulations of a turbulent boundary layer up to Ret = 1410,” 

J. Fluid Mech. 187, PP 61-98. 1 

Stephens, C.A. and Crawford, M.E. 1990 “An investigation into the numerical prediction of ® 

boundary layer transition using the K.Y. Chien turbulence modeling,” NASA CR-185252. s 

Yang, Z. and Shih, T.H. 1991 “A k - e modeling of near wall turbu- lence.” Proc. 4th Int. : - 

Sym. on CFD, UC Davis, ppl305-1310. = 


6 






I!ll nillll IIRRIII1I 111 111 hi 




















REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202 4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0186), Washington. DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 

March 1992 


A k-z Calculation of Transitional Boundary Layers 


6. AUTHOR(S) 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


WU -505-6 2-21 


Z. Yang and T.H. Shih 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-6939 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 105604 
1COMP-92-08; CMOTT-92-05 


11. SUPPLEMENTARY NOTES 

Z. Yang and T.H. Shih, Institute for Computational Mechanics in Propulsion and Center for Modeling of Turbulence 
and Transition, NASA Lewis Research Center (work funded under Space Act Agreement C-99066-G). Space Act 
Monitor, Louis A. Povinelli, (216) 433-5818. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Unclassified -Unlimited 
Subject Category 34 


13. ABSTRACT (Maximum 200 words ) 

A recently proposed k-z model for low Reynolds number turbulent flows (Yang and Shih, 1991) was modified by 
introducing a new damping function / . The modified model is used to calculate the transitional boundary layer over a 
flat plate with different freestream turbulence levels. It is found that the model could mimic the transitional flow. 
However, the predicted transition is found to be sensitive to the initial conditions. 


14. SUBJECT TERMS 


k-t modeling; Boundary layer transition; Turbulence 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 

12 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 






















